Cosmological perturbations in multiple-field inflation 
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We analyze cosmological perturbations to the linear order in the context of inflation with an arbitrary 
number of scalar fields. The fields take values on a non-trivial manifold with a positive definite metric 
and are non-minimally coupled to Einstein gravity. The perturbations are decomposed into three 
different types. The scalar-type perturbations are presented in a gauge-ready form without fixing 
the temporal gauge condition, as well as in terms of gauge-invariant variables. The gauge-ready 
method enables us to impose different gauge conditions which are most suitable to the problem at 
hand. We quantize the scalar perturbations and obtain the solutions. 
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I. INTRODUCTION 

It is well known that the cosmological inflation, an 
epoch of accelarated expansion, provides a causal mecha- 
nism for the generation and evolution of large-scale struc- 
ture formation in the Universe Inflation explains 
a number of puzzles of the Big-Bang theory, such as ho- 
mogeneity, the isotropy of the Cosmic Microwave Back- 
ground Radiation (CMBR) and the flatness of space- 
like sections. As an added bonus, it connects cosmology 
with high-energy physics, thus forming a cosmic labora- 
tory where one may probe physics beyond the Standard 
Model. The very high accuracy CMBR data recently ob- 
tained by the WMAP satellite |4j provide a new impetus 
to compare the predictions of cosmological inflation with 
observations and hopefully to discover new physics in the 
very high energy regime. 

At the simplest level, the inflationary scenario is im- 
plemented by assuming that the matter is described by 
a single scalar field, the inflaton, which is a special case 
of a perfect fluid 0, ■ As the early Universe under- 
goes inflation, quantum fluctuations of the scalar fields 
are generated which become classical after crossing the 
event horizon. During the decelaration phase they re- 
enter the horizon and seed the observed density pertur- 
bations. Cosmological perturbations in a single field in- 
flation has been thoroughly investigated in the past 0,13- 

Despite initial successes, the single component infla- 
tion also has its drawbacks. It was realized quite early 
that, in its original form, the inflationary scenario suf- 
fers from what is called the graceful-exit problem 0,0- 
Linde 7] showed that in order to achieve sufficient infla- 
tion consistent with the observed density perturbations, 
before the Universe exits from the inflationary epoch, 
one requires at least two scalar fields without modifying 
Einstein gravity, and without sacrificing natural initial 
conditions. 



* Electronic address: joy®veccaLerriet^inJ^ 
t Electronic address: 



gautam Qtheory.saha.ernet. in 



Of course, there are other motivations for incorporat- 
ing multiple scalar fields contributing to the dynamics 
of inflation. When constructing models of inflation in- 
spired by particle physics theories such as low energy 
effective supergravity derived from superstrings, one ob- 
tains many scalar fields (see || for a recent review) . Thus 
it is necessary to have a general framework for handling 
cosmological perturbations in a situation where the mat- 
ter sector consists of an arbitrary number of scalar fields. 
A method for treating density perturbations in multi- 
component inflation was proposed in |Ioj . but see also 

[Ulil0 

The study of cosmological perturbations was initiated 
by Lifshitz jlj| in 1946 when he analyzed hydrodynami- 
cal fluid perturbations in Einstein gravity. He assumed a 
particular gauge which is now known as the synchronous 
gauge. This gauge does not completely fix the gauge de- 
grees of freedom and the spurious gauge modes have to 
be properly sorted out in order to obtain correct results. 
Later on, the zero-shear gauge was used by Harrison [l5| 
and the comoving gauge by Nariai [T^j. However, it was 
the seminal paper by Bardeen |T^ which helped put cos- 
mological perturbations on a proper footing. He intro- 
duced a number of gauge-invariant variables to the linear 
order, in terms of which the perturbations became much 
simpler to analyze. Reviews of cosmological perturba- 
tions may be found in and . 

A different approach to cosmological perturbations was 
elaborated b y H wang and colleagues in a series of pa- 
pers 0, |2(J, T2lL l22fl , following a suggestion by Bardeen 
[23 , that rather than imposing a particular gauge condi- 
tion right from the beginning, it is often advantageous to 
express the perturbations without specifying any gauge. 
This then adds the flexibility of adopting different gauge 
conditions at a much later stage, depending upon the na- 
ture of each problem. Moreover, it becomes easy to re- 
late results between various gauge-dependent and gauge- 
invariant techniques. This approach has been termed the 
gauge-ready method. 

In this paper we apply the gauge-ready method to an- 
alyze perturbations in cosmological inflation driven by 
multicomponent scalar fields. This paper is organized 



Typeset by REVTeX 



A The scalar fields 



2 



as follows. In Section [H] we set up the equations de- 
scribing multicomponcnt scalar fields with a non-trivial 
field metric coupled non-minimally to Einstein gravity. 
We then introduce a set of basis vectors, using Gram- 
Schmidt orthonormalization, which enables us to disen- 
tangle multiple-field effects from single-fields ones. The 
background equations, metric perturbations and the per- 
turbed order variables are presented in Section IlIII Here 
we also discuss briefly the issue of gauge transformations 
as applied to cosmological perturbations. In Section HVI 
we introduce the gauge-ready approach to cosmological 
perturbations in the multiple-field inflation scenario. We 
present the perturbation equations in the gauge-ready 
form, as well as in terms of gauge-invariant variables de- 
rived from the gauge-ready equations. Slow-roll variables 
in the context of multicomponent inflation are presented 
in Section We proceed to apply canonical quantiza- 
tion to the density perturbations. The solutions to the 
equations for quantized perturbations are then derived 
to the first order in slow-roll. A brief discussion of vector 
and tensor perturbations is also presented. We conclude 
in Section IVll 

II. PRELIMINARIES 

A. The scalar fields 

In this Section, we explain our notation and set up the 
basic equations needed for our analysis. As our start- 
ing point, we consider Einstein gravity coupled to an ar- 
bitrary number of real scalar fields. We write the La- 
grangean as 

2 = ^U^R-^Q-d^-V^ 



For any two vectors A and B, we define the inner product 
and the norm as 

A B = A f B = A T GB = A a G ab B b , 

\A\ = v/M, (4) 

respectively. Here A* is the cotangent vector such that 
(A^) a = A b Gba- We also introduce the covariant deriva- 
tive V a on M acting upon a vector A as 

V b A a = A\ + T a bc A c , (5) 

while, the covariant derivative on A with respect to the 
spacetime is 

-DuA'sdn + T&d^A . (6) 



9 (J^R - \g^d^ T Gd v 4> - n<£)) • (i) 



Here R is the scalar curvature, Kq = 8ttG, and we set 
c = 1. For the scalar fields we use a vector nota- 
tion, 4> = (0 a ), where the indices a, b, c, ... — 1,2,3, 
. . . , N label the iV-components in field space. Further, 
g = det(g PiL ,), and /i, v, . . . denote the spacetime indices. 
For repeated indices, the summation convention applies. 
The second quantity within the parentheses of Eq. 
represents a nonlinear sigma-model like non-minimal ki- 
netic term. Such a kinetic term appears in various mod- 
els of high-energy physics @- Also V(<p) is an arbitrary 
scalar potential. 

Following the authors of Jnj, we can interpret the 
scalars <j> as coordinates (tfi a ) on a real manifold M. in- 
duced with a symmetric Riemannian metric G having 
components G a b in the field space. The field metric is 
chosen to be positive-definite so that the Hamiltonian 
is bounded from below. The special case of minimally- 
coupled fields corresponds to the situation G a b = S a b- 
From the components G a b we can define the connection 
coefficients T% c in the usual manner, 



^bc ~ -G ad (Gbd,c + G c d^b — Gbc,d) ■ (2) 



The curvature tensor on M is introduced in terms of the 
tangent vectors B, C, D: 



(3) 

I 

It should be noted that the covariant derivative reduces 
to the ordinary derivative when it acts upon a scalar. 

By varying the action Q with respect to g^ and <fr, 
we obtain the gravitational field equation, 

=T^= ■ d v <t>-8» QdV • d x 4> + V^j , (7) 

and the equation of motion for the scalar fields, 

gT (vX-^)dx4>-G- 1 V T V = 0, (8) 

where G^ v and are Einstein and energy-momentum 
tensors. 

It is often convenient in our analysis to represent the 
scalar fields as effective fluid quantities. We conclude 



I 

[R(B, C)D] a ee R\ cd B b C c D d = (r£ d>c - T a bc4 + Tt d T a ce - Tt c T a de ) B b C c D d . 
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this Section by covariantly decomposing the energy- 
momentum tensor into fluid quantities using a time-like 
four- vector m m normalized as u^u^ = — 1: 

T a f3 = fJ-u a Uf3 + ph a p + q a up + q fj u 

H = T a/3 u a u , p = ^T a/3 h af3 , q a = -T^u^hl, 



n a f3 = Tjsh^hp — ph a p. 



(9) 



Here p,, p, q a , and ir a p are the energy density, pres- 
sure, energy flux, and anisotropic pressure, respectively; 
hafi = g a /3 + u a U/3 is a projection tensor based on u a vec- 
tor, q a u a = = 7r Q/ 3, and 7r" = 0. The decomposition 
given above is in the most general form. Indeed, for a 
multicomponent scalar field, we have 



q a = = 7T Q/ 3. 



(10) 



Equations J7J and (JSJ with IjlUI) provide the fundamental 
expressions required for describing cosmological inflation. 



B. Basis vectors 

We continue our discussion of fields on a manifold by 
introducing a set of basis vectors [l^, which will prove 
to be useful in our analysis. First note that an arbitrary 
tangent vector A on M can always be expanded in terms 
of a set of basis vectors {e( a )} as A = A^e^ with 
e (a) ' e (b) — Gab- However, a different set of basis vectors 
generated using Gram-Schmidt orthonormalization turns 
out to be more convenient. 

From the vector <f> we can construct a set of N linearly 
independent vectors 4>( 2 \ ■ ■ ■ , </> ( - JV - ) }, where, 







(1) = 



0, 



<j> (n>2). (11) 



Let ei = tf^ 1 ' /\(fr- l > | be the first unit vector along the 
direction of the field velocity <fi. Define the second unit 
vector e 2 to be along that part of the direction of the 
field accelaration T> t (f> which is normal to e\\ 



e 2 



<ft(2) _ ( Cl . 0(2)) ei 
|0(2)_( ei .0C2)) ei | 



(12) 



It is obvious from Eq. I|12[l that e\ ■ e 2 = by construc- 
tion. A repetitive application of the Gram-Schmidt pro- 
cedure then generates a set of mutually orthonormal vec- 
tors {e„}, which span the same subspace as the vectors 

Introducing the projection operators P n and P^, 
which project on e„ and on the subspace perpendicular 
to {ei, . . . , e n } respectively, we may then write a general 
unit vector e„ as, 



i^_i0 (n) r 



(13) 



where, 



and we define, 



En 



Pi 



eie{, 



p i 



1 



(14) 



(15) 



Note that when the denominator in Eq. (|13l) vanishes, 
the corresponding basis vector does not exist. 

Using the fact that P" + P 1 - = 1, we can decompose 
any vector A in directions parallel and perpendicular to 
the field velocity: 



A = A" + A 1 - = (P 11 +P^)A 
= e x (ei • A) + e 2 (e 2 • A). 



(16) 



For the special case of just one field, e.\ by definition sim- 
ply reduces to the normalized scalar |. Hence, 
from Eq. I|12|) , e 2 vanishes identically, and so do all other 
basis vectors. Thus the decomposition l|ltjfl enables us 
to distinguish between single-field contributions, where 
only ei survives, from multiple-field ones. 



III. THE PERTURBED UNIVERSE 

A. Metric perturbations 

We know that the observed Universe is not per- 
fectly homogeneous and isotropic. Assuming that the 
inhomogeneities are small enough, we can then treat 
the deviations using perturbation theory. In this pa- 
per we consider linear perturbations of the homoge- 
neous and isotropic cosmological space-time described by 
the Friedmann-Robertson- Walker ( FRW ) model. We 
choose the line-element to be 



ds 2 = -a 2 (1 + 2A) dr] 2 - 2a 2 B. l dx i 



+a 2 (<? 



(3) 

ij 



+ 2C ij )dx i dx\ 



(17) 



where a(t) is the scale factor, dt = adrj, and indices 
i,j, . . ., run from 1 to 3 labelling the spatial compo- 
nents. The perturbed order variables A(t,x), Bi(t,x), 

and Cij (t,x) are based on the metric of the 3-surfaces 
of constant curvature K = 0, ±1. These are general func- 
tions of space-time which characterize the linear cosmo- 
logical perturbations. Here t and rj are the comoving 
and conformal times respectively. We denote a deriva- 
tive with respect to comoving time by =9* and one 
with respect to conformal time by / = d n . The Hubble 
parameters in terms of comoving and conformal times are 
defined as H = a/a and Ti. = a' /a = aH. 

As is evident from Eq. I|17|l . the metric is decomposed 
into a background part, plus a perturbation. Correspond- 
ingly we can decompose the scalar field as 



(p(t,x) = <j>(t) + 6<j>(t,x) 



(18) 
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where the perturbation Sep = {5(j) a ) is a tangent vector on 
M , while the energy- momentum tensor is decomposed as 

T° = -n = -(fi + Sfi), 

T°i = - [qi + (fJb + p)ui\ = (fi + p)vi, 



T)= p S)+tt) = {p + Sp)S)+TT 



(3)z 



(19) 



The barred entities denote background variables. For 
notational simplicity we shall ignore the overbars unless 
required. In Eq. i|19|) . is the frame- independent flux 



variable, and 7T j are based on . 

From Eqs. (JZJ) and JHJ, the equations for the back- 
ground can be written as 



,(3) 



(C 



the tensor-type perturbation eff is transverse-traceless 
- Q = Both the vector and tensor per- 

turbed order variables are based on gj? '. We define A 
as a comoving three-space Laplacian, and introduce the 
following combinations of the metric variables, 



X = a(/3 + aj), K = 3(Ha-ip - x ), 

a 2 - 

$0) = +aC {v) . 



(26) 



It is convenient to separate the temporal and spatial 
aspects of the perturbed order variables by expanding 
them in terms of harmonic eigenfunctions Q^ s ' v,t '(k;x) 
of the generalized Helmholtz equation [171 Il8| : 



tj2 1 2 



K 

- (20) 



1 

H -- 
R = Q[2H 2 



K 



1 



4(^ + p) + '- = -^1 M 2 + ^> ( 21 ) 



a 
K 



K 

H+ — 

ar 



V t cp + 3Hcp + G 1 V J V = 0, 
/i + 3ff (fi + p) = 0. 



(22) 

(23) 
(24) 



Taking the Gq and G\ - 3Gq components of Eq. J7J) yield 
Eqs. ((201 and f2*T|) respectively, while Eq. (JHJ) leads to 
Eq. 1)23(1 . Equation (|24|l follows from the conservation 
of the energy-momentum tensor. We shall ignore the 
cosmological constant A in our work; nevertheless it can 
be easily included by making the replacements (J, — > /J, + 
A/kq and p — > p — A//tg- Note that we have explicitly 
retained Jf(= 0, ±1), and only at a later stage shall we 
set K = 0. 



B. Scalar, vector and tensor decompositions 

In order to make further progress, it is customary to 
decompose the perturbed order variables into scalar-, 
vector-, and tensor- type perturbations. To the linear or- 
der, each of these three perturbations decouple from one 
another and evolve independently. Accordingly, the met- 
ric perurbation variables A(t,x), Bi(t,x), and Cij(t,x) 
may be decomposed as 



A = a, B 

- J3) 



= fc + B 



(v) 

1 

(v) 



-.(*) 



(25) 



In this and the following, the superscripts (s), (v) and 
(t) will indicate the scalar-, vector- and tensor- type per- 
turbed order variables. The vertical bar represents a 
covariant derivative with respect to (/>•' and the round 
brackets in the subscript imply symmetrization of the 
indices. The scalar metric perturbations are then given 
by a, (3, 7 and (p. The transverse-type vector pertur- 





bations B\ v) and d v) satisfy B {v)i 



C (V) 1 while 
I* 









Q« = 


^• + - 
k 2 3 




q( v ) \i 






Q (t) \k 




QW = 


Q (t) \j 


= = Q {t) \. 





-s (s) , 



(27) 



Here fc is the wave vector in Fourier space and fc = |fe|. 
We can then write the scalar-type perturbed order vari- 
ables as a(t,x) = a(t,k)Q^(k;x), with similar ex- 
pressions for /?, 7 and cp. The vector- and tensor- 
type perturbations are expanded as £?| = B^q[ v \ 
c (v) _ and C (t) = C (t) g» j n each f these 

harmonic expansions, a summation over the modes of the 
eigenfunctions is implied. In particular, the perturbed 
scalar fields have the expansion 



8<p(t,x) = 5<p(t,k)QW(k;x 



(28) 



From Eq. I|28[) we derive the important conclusion that 
the scalar fields on their own can generate only scalar- 
type perturbations. We also note that to the linear order 
in perturbations, the form of the equations in configu- 
ration space are identical to the corresponding ones in 
Fourier space. Consequently, for maintaining notational 
ease, we do not distinguish between the two. 

In a similar spirit, the fluid variables Uj and 7T can 
be expanded in terms of the harmonics as 

tt (3) ;. = n^Q {s) ; + tt^qW; + w^q% (29) 

while the energy-momentum tensor in Eq. (|19|) has the 
expansion 

T° = -p, = -(p + Su), 

T° i = -l(n.+p)vf + (u. + p)vMQV, 

T) = (p+ 8 P )5) + q (s ) + tt<"> q!^; + tt^ Q W ] . 

(30) 
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For a Universe having the matter sector composed ex- 
clusively of scalar fields, the quantity in Eq. Ijl9|l 
vanishes identically. We then have to the perturbed or- 
der, 



Sfj, = <fi ■ V t ( 

Sp = 4>- v t i 

. a 
{fi+p)v- = 



(fi - a<|0| 2 + VV ■ 5<p, (31) 
- a|0| 2 - VV • 5<f>, (32) 
4> ■ 5(f>, (33) 



where we have written v = for simplcity. It is also 
convenient to decompose Sp into an adiabatic part c 2 (5/i, 
and an entropic perturbation e: 



where c 2 
velocity. 



dp = c 2 S^i + e, (34) 
• p/fi may be interpreted as an effective sound 

C. Gauge transformations 



As mentioned previously, our Universe shows depar- 
tures from ideal homogeneity and isotropy. When the 
deviations are small enough, one can choose a fictitious 
background geometry and consider perturbations about 
it using infinitesimal coordinate transformations. The 
change in correspondence between the background and 
perturbed space-times represented by coordinate shifts 
is called a gauge transformation. Now, since relativistic 
gravity is invariant under coordinate transformations, the 
perturbations are not unique. There are different ways of 
mapping between the background and perturbed parts. 
This leads to what is called gauge degrees of freedom in 
the context of cosmological perturbations. 

Below we briefly summarize the transformation prop- 
erties of various quantities to the linear order in pertur- 
bations HlHIIi. Under a coordinate shift i» = x^+^, 
the scalars, vectors and tensors transform as 



f)r a 

A(£ x ) = A(x x ), A^(x x ) = —A a {x x ), 



(35) 



so that 



A{x x ) = A(x x ) - A, a ? , 

A[i{x ) = A^,{x ) A^^a^, "^er£,/^5 

Afj, v (x ) = Ajj, v {x ) A^y^a^, ^A/j^u^^y (36) 

Writing the temporal part of £ M as £° = a _1 f* and de- 
composing the spatial part as = a" + Q , where 
is based on satisfying £ |* = 0, we find from 
Eq. H36JI that the metric and matter variables transform 



to linear order as: 



a 



a. 



Gauge 


Gauge 
condition 


Temporal 
component 


Synchronous gauge 


ct = 


o 




i 


Zero- shear gauge 


X = 





e= 





Comoving gauge 


v/k 


EE 


e= 





Uniform-curvature gauge 


if EE 





e = 





Uniform-expansion gauge 


K EE 





e = 





Uniform-field gauge 


S(j) ■: 


= 


e = 





Uniform-density gauge 


SfJ, E 


eO 


e = 





Uniform-pressure gauge 


5p E 


5 


e = 






TABLE I: Common gauge conditions 



7 = 7- -6 <p = (p-H?, X = X-?, 
a 



K = k+ \'iH 



A 

7i2 



V = v C 

a 



Sp, = Su.— Sp = Sp — 8<p = 8<j) — (fit; 

~.(s,v,t) _ _(a,i),t) ✓=»(*) _ r {t) 



(37) 



It is immediately obvious from Eq. (|37|l that the tensor- 
type perturbations are gauge- invariant. For the special 
case of scalar-type perturbations to the linear order, fix- 
ing the temporal part £* of the gauge transformation 
leads to different gauge conditions. Table [I] summarizes 
some common temporal gauges. 

From the first equation in (|37Jl . we see that the gauge 
transformation of a involves the term Therefore the 
synchronous gauge condition a = fixes £ upto a con- 
stant of integration, leaving a spatially varying residual 
gauge mode For the remaining gauges in Tabled 

the temporal gauge mode is completely determined. 

We conclude by giving a few examples of gauge- 
invariant combination of variables: 



fx = f - H X, a x = a-x, v- 
Sfi x == Sfi fix, $Px =Sp- px, 



k 

v X, 

a 



84> - jjtp - 



S(fi x =S(fi- (fix, 8(fi v 

aH a . 

LPv = Lp k~ Vl ^" = d M - r/w- 



H 



9d<t>, 



(38) 



Thus, in the zero-shear gauge, also known as the lon- 
gitudinal, or conformal Newtonian gauge, x = is the 
gauge condition. We then have from Eq. (|38p. ip x = <p, 
a x ee a, and S(fi x = S(fi. Similarly, in the uniform- 
curvature gauge, it follows that 8(fi v = S(fi which in turn 
is equivalent to — (0/ H)ips4> in the uniform-field gauge. 
In the notation of Q , our ct x and ip x correspond to their 
$ and — W respectively. 
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IV. SCALAR PERTURBATIONS IN 
MULTIPLE-FIELD INFLATION 

A. Perturbation equations in the gauge-ready form 



We shall now briefly discuss the gauge-ready approach 
introduced in 0, |2£j, |2ll ■ As is well known in the 
theory of cosmological perturbations, a judicious choice 
of gauge conditions often simplifies the mathematical 
structure of a particular problem. For example, density 
perturbations with hydrodynamical fluids are most con- 
veniently treated using the comoving gauge, while the 
uniform-curvature gauge simplifies the analysis of per- 
turbations due to minimally coupled scalar fields. Since, 
in general, we do not know the optimal gauge condition 
beforehand, it becomes advantageous to express the per- 
turbations without imposing a specific temporal gauge 
condition. In other words, we write the governing equa- 
tions in the gauge-ready form, which would give us the 
freedom to choose different gauge conditions, as adapted 
to the problem, at a later stage in the calculations. Once 
the temporal gauge mode is completely fixed so that no 
further gauge degrees of freedom are left, the resulting 
variables would then be gauge- invariant. Moreover, when 
a solution in a particular gauge is known, we can then 
easily derive the corresponding solution in other gauges, 
as well as in gauge-invariant forms. This is the basic 
concept of the gauge-ready method. 

To implement this gauge-ready strategy, it is most con- 
venient to derive the perturbed set of equations from the 
(3+1) ADM [Hi, and the (1+3) covariant [H formula- 
tions of Einstein gravity. A complete set of these equa- 
tions may be found in the Appendix of Ref. |l9j| . In this 
Section we write the equations for scalar-type perturba- 
tions in the gauge-ready form. 
Definition of k: 



1 Ik 2 
<p = Ha - -k + x^X- 
3 3 a 2 - 



(39) 



ADM energy constraint [G% component of the field equa- 
tion): 



k 2 -3K 1 a . 
— 2 <p + Hk= --k 5(j,. 

ADM momentum constraint (G? component): 

k 2 -3K 3 ,. N a 
K a 2 X= 2 K °^ + P >k V ' 

ADM propagation(G* — |<KG| component): 

a 2 

X + Hx~ a - ip = kI—tt {s) . 
Raychaudhuri equation (G\ — Gq component): 

k + 2Hk + [3H - — a = -Ko(5(i + 3Sp) 



(40) 



(41) 



(42) 



(43) 



Equation of motion for scalar fields: 



A 



A 



Vf + 3HV t - + M z j S<j> = ^a~3<p-^x)(t> 

- 2aG- 1 V T V. (44) 



Energy conservation: 

Sfi + 3H(Sfi + dp) = ((j, + p)(k- 3Ha - . (45) 

Momentum conservation: 
[a 4 (^ + p)v] 



a 4 + p) 



1 



M + P 



. 2fc 2 
Sp ~3^ 



3K 



(46) 



In the above equations, Sfj, and Sp are given by 
Eqs. (jBTjl and (j^2*|l respectively, while 



M 2 = G 1 V T W - R{<t), <p). 



(47) 



Note that these equations are valid for any K, and for a 
scalar field, = 0. 

Equations (|39|l - l|46|l . together with the background 
equations l|2U|) - (|24|l . and the perturbed order variables 
for the scalar fields H31|) - (|33|) . provide a complete set of 
equations for analyzing scalar-type cosmological pertur- 
bations with multicomponent scalar fields. As we have 
not chosen a specific gauge so far, Eqs. (|39|I - H46|) are 
therefore in the gauge-ready form. This allows us to 
impose any one of the available temporal gauge condi- 
tions, which would then fix the temporal gauge mode 
completely, leading to gauge-invariant variables. 



B. Perturbation equations using gauge-invariant 
variables 



In order to illustrate the gauge-ready method, we 
derive some useful expressions in terms of the gauge- 
invariant variables introduced in Section [ill CI 

From Eqs. flUI) and lUTJ) we obtain 



k 2 -3K 



-fx 



-Kq8{i v . 



Eq. 142|) can be written as 



fx = - K o^ 



Eqs. and lead to 



fx ~ Ha x : 
Eqs. gSl,® with @TJ yield 



1 2/ \ u 

-K {n + p)-v x , 



(48) 



(49) 



(50) 



8fi v + 3H6[i v 
k 2 - 3K 



(51) 
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while Eqs. (02) and give 
v x + Hv x 

Sp v 



2lf_ 
3 



3K vr( 



/i + p 6 a* fi- 
Combining Eqs. (|48H - 152I) we can derive 



(52) 



</3 x + (4 + 3c 2 s )H^ 



( M cj-p)-2(l+3cfW 



= _I K 2 C e _H (.) 
2 I 3 



1 

'2 Ko_ H" 



2ff 2 



- p k 2 



(53) 



where we used Eq. <|34|) . For the explicit forms of the 
gauge-invariant variables used in these equations, see 
Eq. (|38|l . From Eq. (|49|l we can draw the important con- 
clusion that, for scalar-fields, a x — —(p x , since ir^ = 0. 
Using this result the equation of motion for scalar fields 
becomes 



(V 2 + 2HV n - A + a 2 M 2 ) d(f> x = -4sp' x <f>' 
+2a 2 f x G- 1 V T V. 



(54) 



It is also convenient to re-write Eqs. I|5U|I and l|53|) for the 

case of scalar fields as 



(55) 



ip' x + Hip x 

tp x + 6H<p' x -A<p x + 2 [H' + 2(H 2 - K)] <p x 

= n 2 a 2 VV ■ 5cf> x , (56) 

where we used Eq. (|33|) . and the relations 

e = Sp- c 2 Jn = Sp x - c 2 Jn x , 

(1 - 4)5u. x -e = 8fx x - 5 Px = W ■ 8<t> x . (57) 

Eq. (|55|) is often called the constraint equation. Note that 
we have written Eqs. I|54|) - (|56|l in conformal time. These 
equations contain most of the physics related to infla- 
tionary cosmological perturbations. They are expressed 
in terms of gauge-invariant forms of the variables, and 
from the discussion at the end of Section IIII CI we see 
that they retain the same algebraic forms in the zero- 
shear gauge. 

We end this Section by expressing Eq. (|5fc>|) in a differ- 
ent way. Observe that according to Eq. (|16[) . 8<p x may 
be decomposed into components parallel and perpendic- 
ular to the field velocity, 8(p x — 8<p x + 8<p x . Using the 
background equation l|23|l . the constraint equation 1)55(1 . 
and the fact that |</>'|'|</>'| = (2? v 4>') ' 4>\ we can write 
Eq. ® as 



2 [H- 



n ~ n W\ 

4(V v <f>') ■ 8cj> A 



- 2K 



(58) 



Following our discussion in Section III Bl we know that 
the perpendicular component of field perturbation van- 
ishes when there is only one field. In this case, the right 
hand side of Eq. (|55j) vanishes, and the resulting equa- 
tion is well known in the theory of single field inflationary 
perturbations 



V. SOLUTIONS OF THE PERTURBATION 
EQUATIONS 

A. Slow-roll variables 

To proceed further with our analysis, we make the as- 
sumption that the Universe has undergone inflation to 
complete flatness, so that henceforth we can set K = 0. 

We introduce the functions, 



H 0< 2) 



(59) 



known as the slow-roll variables. Using Eq. I|16[l . V ^ s 
decomposed into parallel and perpendicular components: 



?7" = ei • t] = 



V t <P-(f> 



V = e 2 • V = 



H\<f>\ 2 ' ' - ' H\<j>\ 
The standard slow-roll assumptions are 

v 1 - = o(0, 



(60) 



0(C), 17"= 0(C) 



(61) 

for some small parameter £, with e, \fer^ and y/en 1 - much 
smaller than unity. If in an expansion in slow-roll vari- 
ables we neglect terms of order 0(( 2 ), we claim that 
expansion to be of first order in slow-roll. Thus terms 
with e 2 , ery", etc. are of second order. Note that the 
definitions (|59|) remain valid irrespective of the slow-roll 
assumptions. 

We present some useful relations involving the slow-roll 
variables: 



H' = H 2 (l-e), nLL = H(l + T)"), 
V v <f>' = H\<f>'\(ri + ex) = n^ 1 V2H 2 V~<r) + ex), 



n 2 e= -K 2 \ct>' 



l\2 



e' = 2We(e + r?H). 



(62) 



B. Analysis using gauge-invariant variables 

In order to solve the system of perturbation equations 
l)54p. (|55|l and (|58(l . we shall find it convenient to intro- 
duce the variables 

q = a (s$ x - —cp x j =a(64>- —<p) , (63) 
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where the second equality for q in Eq. (|63(l follows from 
Eq. 1138(1 . Indeed, q is gauge- invariant, and is a natural 
generalization of the single field Sasaki-Mukhanov vari- 
able 0. 

Using the slow-roll variables 1(59(1 . together with some 
of the relations JfJ^J), the constraint equation l|55|l can be 
written in terms of q as 



1 2,1 1 

2 a 



^+tt(l + eV x = -~^0'-± (65) 

From Eq. (|54() . the scalar field perturbations satisfy 
V\q- (A- H 2 n)q = 0, (66) 



where 



n = - (2 - e)l - 2e ((3 + e )P« + erf + ryej 

(67) 

and use has also been made of Eq. H65|) . The correspond- 
ing Lagrangean Jz? follows from Eq. 1(66(1 : 



S 



X - J (VrfV v q + q^A~H 2 n)q) ^)dnd 3 x. 



(68) 



Here is the determinant of the metric g\j of the 3- 
surfaces of constant curvature K = 0, see below Eq. (117(1 . 

The equation of motion for u is obtained by substitut- 
ing its definition into Eq. (|5*5j) : 



Am - 

n 



-u = K,lnn ± q 2 , q 2 = e 2 - q, 
k 1 



a|0'| V2ay/I' 



(69) 



For later use, we also express Eq. (I^Sll in terms of u 
and q as 



(W , 

1/0 



1 



-gi, q 1 = e 1 -q. 



(70) 



Differentiating Eq. I|7U|) once with respect to the confor- 
mal time and using Eq. I|69|l . we obtain the relation 



2 I «i 



(W 

1/0 



<7i - K 7i?7 92 = Au. 



(71) 



Although the equations (|fJ5|) . (|rJrj|) and l(fi9"|) have been 
expressed in terms of the slow-roll variables, they are ex- 
act, and no slow-roll approximation has yet been made. 
Observe that, to the leading order in slow-roll, the per- 
turbation variables q and u decouple, whereas at first 
order, mixing between these occur. 



C. Quantization of the scalar perturbations 

We now study the quantization of the density pertur- 
bations described by the Lagrangean in Eq. JBSJl. We 
start by introducing the matrix Z mn defined as 

rt 

where the second equality follows from D,(e m • e„) = 0. 
Thus Z is antisymmetric and traceless (TrZ = 0). Upon 
expanding q = q rn e m , using the basis {e m }, it follows 
from Eq. JBBJl, 



i (<?' + UZq) T {q' + HZq) + ^q T (A - H 2 n)q, (73) 



where (Q) mn — ej n tte n , and for notational ease, we have 
suppresed the indices to, n. It will prove convenient to 
reduce the Lagrangean 1)73(1 to the canonical form. We 
redefine q using a new matrix R as 

q{ri) = R{rj)Q{rj), R' + HZR = 0, f2 = R T flR. (74) 

From the equation of motion ((74(1 for R, it follows that 
R T R and deti? are constants, so that R represents a rota- 
tion. Without any loss of generality, the initial value of R 
may be chosen as i?(^o) = 1- Substituting the variables 
defined in Eq. (J73J) into Eq. (U3 yields 



% = t:Q Q' + t:<2 (A - n 2 n)Q. 



(75) 



To proceed with the quantization, we employ the 
canonical quantization procedure to the Lagrangean 17511 . 
The momentum II canonically conjugate to Q is 



U(tj,x) = dJ?/0Q /T = Q'( V , x). 
The Hamiltonian is then given by 

jr = in T n - \q t (a - n 2 n)Q. 



(76) 



(77) 



The canonically conjugate variables (Q, II) are promoted 
to quantum operators (Q, TV) satisfying the commutation 
relations 



[a T Il(r,,x),0U(r ] ,x')}=O 
[a"Q(n, x),f3fl(rj, x')] = ia T {3S(x - x'), 

where the delta function is normalized as 

x')y/g~Wd 3 x, 



[a T Q(n 7 x),PQ(7 h x')} 



5{x 



(78) 



(79) 



and we have introduced the vectors a, f3 with compo- 
nents a m , Pm in the basis {e m } to avoid writing the 
indices to, n in the commutators. Since we are consider- 
ing spatially flat hypersurfaces (K — 0), the operator Q 
may be expanded in a plane wave basis as 



Q = 



d 3 k 

(2tt) 3 / 2 



Q* k {rj)a k e 



ik-. 



Male 



t „ — ik-x 



(80) 
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with a similar expansion for II. It immediately follows 
from Eq. 174|) that q must now be interpreted as the 
operator q with modes 

q k ( v ) = R(v)Q k {v), (81) 

satisfying a mode expansion identical to Eq. JHOJ. The 
creation and annihilitation operators a k and ctk satisfy 

[a T dfe,/?a fc /] = [a T a k ,f3a k ,] = 0, 

[a T a k , /3a{,} = a T (35{k - k'). (82) 

In order that the commutation relations Ij78(l and (|82|1 be 

consistent, the following Wronskian condition must be 
satisfied, 

W{Q k , Ql) = Q' k (v)QUv) - Qk(v)Qk(v) = i. (83) 

From the mode expansion i|80|) and the Hamiltonian (|77|l , 
it follows that the equation of motion for Qk is 



Ql + (k 2 + H 2 VL)Q k = 0. 



(84) 



It may be easily verified using Eq. (jHH that the Wron- 
skian satisfies dW{Qk, Q%}/d-q = 0. 

We also interpret the variable u introduced in Eq. H64|) 
as an operator u, and after performing a mode expansion 
identical to that of Q in Eq. ljSU|) , it follows from Eq. 
that the modes u k satisfy 

u k = Kf ) Hr) J -q 2 k, q%k = {e-2 • e m )q k , 

(85) 

or, equivalently, from Eq. I|71|l . 

K&ir^qik - ^ (q[ k - ^J^-^ik) = k 2 u k , 



qik = (ei • e m )q k . 



(86) 



D. First order solution 



As a prelude to presenting the solution of the pertur- 
bation equations to the first order in slow-roll, we briefly 
discuss the issue of the ambiguity in the choice of the 
vacuum state whe n q uantizing fields in an expanding 
FRW background |2g. In ordinary Minkowski space- 
time, there exists a unique time direction, as well as 
distinct time-invariant positive- and negative-frequency 
modes. However, when quantizing in a curved back- 
ground, there is neither a distinct time direction, nor a 
notion of time-invariant mode. Consequently, there is no 
unique vacuum state either. Hence if we have a positive- 
frequency mode Q^ at times 77 < 770 < 0, with the initial 
vacuum state afc|0) = 0, then at later times r\ > \t)q\, the 
modes will in general be described by a linear superposi- 
tion of positive- and negative- frequency modes Q^ , Q^ , 
related by means of a Bogoliubov transformation: 

Qk(v) = x k{v)Qt( 1 lo) + Hk(v)Qk(vo), 
\\k(v)\ 2 - \Mv)\ 2 = 1- (87) 



For the coefficients Ak and (ik, one often makes the choice 
of the initial values as 



Afe(»7o) = 1, HkiVo) = 0, 



(88) 



for the adiabatic vacuum (or the Bunch-Davies vacuum 
in de-Sitter space), which corresponds to the positive- 
frequency solution in the Minkowski space. 

To proceed further, it is convenient to introduce the 
time rj-H when the mode with wave number k crosses the 
Hubble radius during inflation, so that the relation 



H 



(89) 



is satisfied for each k. Consequently, the inflationary 
epoch can be separated into three regions: the sub- 
horizon region (7i <C k), the transition region (7i <~ k), 
and the super-horizon region (TL S> k). We now discuss 
each of these in turn. 

In the sub-horizon region, we solve Eq. Q84J1 with the 
7i 2 Q, term subdominant compared to k 2 . The solution is 
obtained in the limit k/TL — > 00 for fixed k as 



Qk(ri) = 



I 



Jk(v-Vo) 



2k 



R(vo) 



1. 



(90) 



Since one is usually interested in calculating quantities 
at the end of inflation, this region is therefore irrelevant. 

We consider next the transition region. It will prove 
useful to introduce the time r\- when the sub-horizon 
epoch ends and the transition region begins. In a suffi- 
ciently small interval around 77-7^ we can then apply slow- 
roll to the Eq. I|84|) keeeping all the terms, but taking the 
slow-roll functions to be constant to the first order. The 
initial conditions are chosen as 



Qkbi-) 



(91) 

Integrating the relation for 7i' in Eq. (|62|l with the initial 
conditions for the transition region, we obtain 



(1 - e H )r,' 



{1-enW 



(92) 



so that Tt(rm) = k. Differentiating 9 in Eq. 169(1 yields 
9'/8 = —H(l + e + 77II), which can be integrated with the 
result 



(z) = n [ — 



(1 + 26-K+7)" ) 



, z = kr], 



Oh = -j= ^ , z H = kri H (93) 

The differential equation (|74|) for the rotation matrix R 
can be solved with the initial conditions l|91|l leading to 



R(z) = — 



-(l-e n r Z r 



z_ = fc?7_. (94) 
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Since the time-dependent terms in the matrix Q in 
Eq. (|67f) are of first order, we can take fl = Q(t]t-c) = 
in the transition region. Then the matrix Q is given by 



n = R- 1 {z)n H R{z) = n H -[n Hl z n ]\n— 

Z- 

( z 3 

= Q-h + 3[5n, Zn) In h -In e H 

V z n 4 



(95) 



with 



-3 I 21 







= el- 



- + 2e(ei • e m )(ei • e„) T , 



<5h = S(rm). 



(96) 



Here the second equality is valid to the first order in slow- 
roll, and we have used the notation M 2 = eJ n M 2 e n . 
We also made the assumption that those components 
of a 2 M 2 /7i 2 which cannot be expressed in terms of the 
slow-roll variables are of first order. Because dn and Z-j-c 
are both of first order, we can take (l = in Eq. I|95|) 
to be a first order quantity. 

In order to write the equation for the mode Qk in 
the transition region, we will find it convenient to define 
Qk = RnQk{z) and Cl = Rn&Ry 1 , with Rn = R(zu)- 
From Eq. (|94[l. we have to the first order, Qk{z) — Qk(z), 
while from Eq. I|95l) we conclude that Cl = VL-h within 
a small region around zn- Using the above results in 
Eq. (|84|) . the mode equation for Qk may be written in 
terms of Qk as 



Qk, zz 



1 - 



Q k = 0, v 2 H = -t + 3S H . (97) 



This equation is similar to the one obtained for the single- 
field inflation, except that this is a matrix equation. The 
solution is then given in terms of the Hankel functions of 
matrix valued order vy^, 



Qk(z) = ^[ Cl (k)H^(z) + c 2 (k)H^(z)], 



(98) 



We wish to match the solution in Eq. (|98(l so that in 
the limit k/TL — > oo, the modes approach plane waves, 
Q k (z) = e lz /V2k, see Pty. For \z\ > 1, the Hankel 
functions have the asymptotic forms, 



H^(z) ~ y r 2/{^z)e^-^ +1 / 2 ^' 2 \ 
H™(z) ~ ^2/(^7)e-^ z -^ +1 / 2 ^l 2 \ (99) 



We set ci(fc) = v /^7(4fc)e l(l/H+1/2),r/2 , and c 2 (k) = 0. 
The phase factor of c\ (k) is chosen in order to match with 
Eq. (|90|l at short scales, while the factor of WTtJ (4k) en- 
sures conformity with the Wronskian in Eq. (|83|l . There- 
fore the final solution with the appropriate normalization 



Q k (z) = VWWy {un+1/2)7r/2 ^H^(z). (100) 
It is worth mentioning that that the matrix valued Han- 
kel functions are to be interpreted as series expansions, 
just like the usual Hankel functions. 

We finally discuss the solution in the super-horizon re- 
gion. On super-horizon scales we have \z\ <C 1, for which 
the asymptotic form of the Hankel function is 



H$(z) 



/2^ e - il /2 2 ^-3/2£W^- l 



r(3/2) 



(101) 



so that the asymptotic solution for Qk(z) in the super- 
horizon region is given by 



Qk(z) ~ (l/V2k)e l{u ' H - 1/2)7r/2 2^- 3 ^ 2 



T(v H ) h. 

Z * 



r(3/2) 

-(l/Vuy^-^+^WEH^/zH)- 1 - 5 ", 

(102) 



where 



E n = (l- e H )t + (2- 7 B - In 2)5- 



(103) 



and je ~ 0.5772 is the Euler constant. 

In this region since k/TL — > 0, we can also solve Eq. (|85|) 
ignoring the k 2 dependent term, leading to 



Uk(ij) =u P k + C k 9 + D k 9 



dr)' 
6 2 (r)') 



up k = 



" drf_ 

„ ~o 2 



dn"dK,lHir L q2k, 



(104) 



where Ck and Dk are constants of integration, and up k 
is a particular solution. Note that since 9 is a rapidly 
decaying function, we can ignore Ck compared to Dk . In 
the same approximation, the solution of Eq. (|86() is 

qik = d k (l/9) + 2(1/9) f ' d V '9K 2 Hr ] ± q 2 k. (105) 

J VH 

From Eq. I|70[l we see that the integration constants Dk 
and dk are related by Dk = ^dk- Considering the re- 
gion where r\ is sufficiently close to 77^, the integral in 
Eq. I|105|) may then be neglected, so that using Eq. 193|l . 
we can write q\k = %D f.(l / 0?i) ( z / z-h)^ 1 ■ Taking into ac- 
count the asymptotic solution (|102fl . and the fact that 
q k = (ei • e m ) T qik, we finally obtain, 



D k = -(l/2V2fc)e 



J(uh-1/2+2Sh)tt/2 



9n(ei ■ e m ) T E H - 



(106) 

Thus the integration constant in Eq. (1 104(1 is completely 
determined to first order in slow-roll. Inserting the result 
(|106fl for Dk in (|104fl . and using the relation anH-j-c = k, 
we finally arrive at 
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Uk = ~(2^ ei( ' H ~ 1/2+25 ' HW2 ^ W*«.*)( e i ' e ™) T + B (*«'*)] ( 107 ) 



where we ignored Cfc, and 



A(t H ,t) = -^=£ *' a (^) ' B ^ n ^ = a^fej t dt'a(J^ U(t n ,t), 

U(t H ,t) = 2kI f dt'Hr ] ± x f^^(e 2 ■ e m f R 
Jtu yea 



(108) 



Here Qkn is the value of the asymptotic solution I|1U2|) 
for Qk evaluated at r\ = rju- Observe that the solution 
(I1Q7|I for uu is expressed entirely in terms of background 
quantities and comoving time. This concludes our dis- 
cussion of scalar perturbations in multiple field inflation. 

E. Vector and tensor perturbations 



where and Bk are integration constants. Observe 
that VI/W, v( v \ n^ v ' and appearing in these equa- 
tions are gauge-invariant, see Eq. IpFTjl . It is interesting to 
note that we can derive the equations for vector and ten- 
sor perturbations without taking into account the scalar 
fields. Therefore the presence of scalar fields do not for- 
mally affect these perturbations. 



For the sake of completeness, we now present a 
brief discussion of vector- and tensor-type perturba- 
tions. From the G\ component of Eq. JJJ, together with 
Eq. Jj23J|, we have 



while the condition Tff = yields 



-kir {v) . 



(109) 



(110) 



Equations 1)1091) and (|110fl describe the vector-type, or 
rotational perturbations. Since vector sources are absent 
when the matter sector is composed entirely of scalar 
fields, the vector-type perturbations are therefore irrele- 
vant in the inflationary scenario. 

The equation for the tensor-type, or gravitational wave 
perturbations follows from the G l j component of (7J|: 

CW" + + k 2 C^ = K 2 a\^ . (Ill) 

For scalar fields we have = 0. We can recast 

Eq. iJTTTJl as 

v' t '+ (k 2 -^jv t =0, v t = aC^\ (112) 

which is of a form similar to Eq. Q85|). The solution in 
the large-scale limit is then obtained by ignoring the k 2 
dependent term: 



r i 

A k + B k I —dr], 



(113) 



VI. CONCLUSION 

In this paper we presented a general framework for 
analyzing linear cosmological perturbations in the multi- 
component inflation scenario using the gauge-ready ap- 
proach. Our model consists of multiple scalar fields in- 
duced with a positive-definite, but a general field metric, 
coupled non-minimally to Einstein gravity. The space- 
time metric is chosen as the perturbed FRW world model. 
We gave the complete set of perturbation equations in 
the gauge-ready form, and derived a set of equations for 
gauge- invariant perturbed order variables. We wrote the 
equations governing scalar perturbations using general- 
ized forms of slow-roll variables. We then applied canoni- 
cal quantization to the scalar perturbations and obtained 
the solutions to the first order in slow-roll. We found that 
the asymptotic solutions in the super-horizon region are 
given in terms of Hankel functions. This is similar to the 
single-field inflation case, except that the order of the 
Hankel functions is matrix valued. 

There are a number of possible extensions to our work. 
First, the immediate next step would be to calculate the 
power-spectra and spectral indices in realistic models of 
multiple-field inflation. Second, it would be interesting 
to interface our approach with the CMBFAST [27| or 
the CAMB [28| computer codes and compare with the 
WMAP results. Third, the methods in this paper may 
be extended to include generalized gravity as well as hy- 
drodynamical fluids. Fourth, a systematic investigation 
of adiabatic and isocurvature perturbations in the con- 
text of multicomponcnt inflation can be made using the 
gauge-ready approach. These, and other issues, will be 
presented in a future work. 
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